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GALAH SURVEY: CHEMICALLY TAGGING THE THICK DISK 
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Abstract. The GALAH survey targets one million stars in the southern 
hemisphere down to a limiting magnitude of V = 14 at the Anglo- 
Australian Telescope. The project aims to measure up to 30 elemental 
abundances and radial velocities (~ 1 km s -1 accuracy) for each star 
at a resolution of R = 28,000. These elements fall into 8 independent 
groups (e.g. a, Fe peak, r-process). For all stars, Gaia will provide 
distances to 1% and transverse velocities to 1 km s^ 1 or better, giving 
us a 14D set of parameters for each star, i.e. 6D phase space and 8D 
abundance space. There are many scientific applications but here we 
focus on the prospect of chemically tagging the thick disk and making a 
direct measurement of how stellar migration evolves with cosmic time. 


1 Introduction 


The GALAH survey (2014-2019; see Fig. |T| is the latest major undertaking at the 


Anglo-Australian Telescope (De Silva et al. 2015) It uses the $12M HERMES 


instrument fed by 400 fibres that can be positioned robotically at the prime focus 
(Sheinis et al. 20 14)| The main science goal of the project is to chemically ‘trace’ 
different Galactic components aided by the stellar kinematics with a particular 
emphasis on old stars, i.e. stars that were born (half of the stellar mass) before 
z ~ 1. In each of these components, we propose to explore chemical tagging, in¬ 
troduced as a tool for reconstructing information that has been lost over billions 
years as the Galaxy evolves (Freeman & Bland-Hawthorn 2002)]' The central idea 
is that essentially all stars have been born in homogeneous gas clouds that have 


long since dispersed (Bland-Hawthorn et al. 2010a) Only a small fraction remain 
bound today in the form of star clusters (e.g. open clusters). Thus, in certain 
instances, a star’s siblings may be identified by a unique chemical signature in a 
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Fig. 1. Potential GAL AH survey fields (4300 in total) across the southern sky are shown 
as open grey circles; the magenta circles indicate the fields observed up to November 
2014. The full survey will target about three quarters of the potential fields. 


high dimensional chemical space (C— space). In principle, we can hope to recon¬ 
struct where some stellar families were born, and how they have become dispersed 
throughout the Galaxy. 


Arguably, the most interesting target is the thick disk believed to be more 
than 10 Gyr old. Of all components, we predict that the thick disk will carry the 
strongest signature of chemical tagging. The thick disk has a very well defined 
enhancement in a/Fe at all metallicities (Bensby et al. 2014) This is important 


because it means the stars were likely born in high pressure, bursty star-forming 
regions (Aalto et al. 1995) Massive star clusters have been observed in high- 
pressure turbulent disks at high redshift although their association with the thick 
disk is uncertain (Lehnert et al. 2014) Direct evidence of enhanced a/Fe ratios 
under starburst conditions comes from the small fraction of metals detected in 


x-ray winds from these environs (Martin et al. 2002) 


This has two consequences for chemical tagging: the initial cluster mass func¬ 
tion (ICMF) is flatter in starburst environs meaning that there are fewer low mass 
clusters, and the extreme clusters tend to be more massive than found in more qui¬ 
escent environments. Both of these conspire to make the tagging signal stronger 
in the thick disk, at least in our models. The quiescent outer regions of the thin 


disk are characterised by in situ formation of smaller star clusters (Larsen 2009) 


We now investigate the expected number and size of clusters in the abundance 
space C for different Galactic components in an observational survey like GALAH. 
Specifically, we investigate as to how the ICMF slope 7 and the maximum cluster 
mass m max affects our ability to detect clusters in abundance space. 
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2 Method 


The Galaxy is modelled as consisting of the thin dislf] (2.15x 10 10 Mg), the thick 
disk (3.91 xlO 9 M 0 ), the stellar halo (7.6x10 s Mg) and the bulge. Depending 
upon the geometry and the magnitude limits of the survey, different surveys will 
sample different fractions of the galaxy. Let /sample be the fraction of stars of a 
Galactic population that is randomly sampled by the survey. Let M pop be the 
actual stellar mass of the population, to* is the mean mass of stars for a given 
IMF (here assumed to be 0.42) and / m i x is the fraction of mass of the population 
that can contribute stars to the survey volume. Then 


/sample 


/VP°P 

survey 


/mix-Mpop/m* 


( 2 . 1 ) 


The factor f m ; x is the fraction of mass of the population that can contribute stars 
to the survey volume. If stars are uniformly mixed over the whole galaxy then 
/mix = I- More realistically, we assume stars born in an annulus of width A R 
around the Sun such that 


/mix 


C-AR/2^ R )^ RdR 


Mr, 


( 2 . 2 ) 


In practice, A R is never zero because orbit families develop radial excursions 
during their lifetime. To accommodate radial excursions, we assume / m i x lies in 
the range 0.25 to 1. For GALAH we have N SUIvey = 10 6 stars. We simulate the 


selection function using the Galaxia code (Sharma et al. 2011) and find that 24% 


are thick disk stars and 75% are from the thin disk; 
pc, half are within 1 kpc and 85% within 2.5 kpc. 


16% of stars are within 500 


3 The initial cluster mass function 


The stars are assumed to be born in clusters and their ICMF is modeled as a 
power law (Elmegreen & Efremov 1997), The cumulative distribution is given by 


/(^ i X m i n , 


o'*! - \ 

•^max ^ 1 

t 1+ 7 _ r 1+ 7 / 

X max X min j 


(3.1) 


with —2.5 < 7 < —1, and size x is in the range Xmin to rr m ax- The mean cluster 
size is then 


x = 


x£(x)dx 


fi+7^ ( *^max - Zmtn \ 

\2 + 7/ - x^ J 


(3.2) 


1 The low disk estimate is from the Besangon model which we use for consistency here; our 
new models presented elsewhere have updated most of the Besangon parameters. 
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The full cumulative distribution of the number of clusters above a certain size x 
is given by 


N(> x\^, 

*£min? ^max) 


/ - x 1+1 A M f mix/m* 

y Imax — J ®(T) x min, r7l max ) 


,.1+7 _ 7.I+7 \ 

m a v 1 


1+7 
2 -max % 


1+7 



+1+7 _ 

- 0, Mfmix/m * 

^X ma l - x min 7 ) 


4 Predictions for cluster size distribution 

Given M pop , /’mix» 7 ) £ m in, imax and N smveyi the number of clusters as function of 
size n in survey is given by 7V(> n// samp i e | 7 , £ m i n ,+’ max ). As a concrete example, 
we predict the number of groups in thick disk that can be seen by GALAH. Here 
ATP u ° r P ey = 2.3 x 10 5 stars, M pop = 3.9 x 10 9 M 0 . 

In Fig. [2j we show the cumulative distribution of clusters as a function of their 
size in the survey. Each panel shows results for four different values of 7 . The 
panels differ in the values of / m j x , N s urV ey and ^max- It can be seen that if / m i x 
is large we get fewer big clusters (top two panels). If x max is small, we expect 
fewer big clusters. If 7 is small we expect to see a larger number of big clusters. 
The bottom right panel shows that if the number of survey stars is decreased by 
a factor of 10 , this dramatically reduces the chances of detecting groups with size 
greater than < 10 . 

In Kig. J3j we show the case for all stars. Here 7V s P° r P ey = 10 6 stars, M pop = 
2.55 x 10 1 < 7 Mq. The total number of clusters increase but the maximum size of 
clusters is smaller. This is because as compared to the thick disk case the ratio 
A?u°rvey/Afpop is smaller. 


5 Survey simulations 


The primary motivation for the GALAH survey is the chemical tagging experi¬ 
ment described in § [T] Our goal is to identify debris of disrupted clusters and 
dwarf galaxies. We assume here that all of the disrupted objects whose orbits 
pass through a ±1 kpc-wide annulus around the Galaxy at the solar circle are 
represented within the observable horizon. Simulations show that a random sam¬ 
ple of a million stars with V < 14 will allow detection of about 20 thick disk 
dwarfs from each of about 3000 star formation sites, and about 10 thin disk dwarfs 
from each of about 30,000 star formation sites (Bland-Hawthorn & Freeman 2004) 
Even a few rare outliers will be enormously valuable for tracking stellar migration 
over cosmic time. The specific example of the Solar Family is discussed else¬ 


(Liu et al. 2015) 


where (Bland-Hawthorn et al. 2010a) and the first searches have already started 


The simulated numbers above depend on the details of the initial cluster mass 
function of the disrupted objects and its mass range. Let us assume the ICMF to 
be a power law specified by minimum mass x m j n , maximum mass x max and power 
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Fig. 2. Cumulative distribution of number of clusters as a function of cluster size for a 
simulated thick disk. 






Fig. 3. Cumulative distribution of number of clusters as a function of cluster size for all 
stars, i.e. both thin and thick disk stars. 
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law index 7 . Our expectation is that the thin disc has a steep ICMF (large 7 ) 
and low x max typical of a quiescent star formation history; the thick disk will have 
a shallower ICMF (small 7 ) and large x max representative of the turbulent high- 
pressure discs seen at high redshift. The size of the GAL AH survey is selected to 
probe the slope 7 and x max of the ICMF. 


We simulate the GALAH survey based on the Galaxia code (Sharma et al. 2011) 


and then applied a simple analytical prescription for generating stellar clusters. 
Fig. [4] shows the number of unique clusters identified in this simulated survey as a 
function of initial cluster mass. Setting the minimum requirement of ten member 
stars for a reliable identification of a cluster, we calculate the threshold cluster 
mass (i.e. the lowest-mass cluster we expect to recover ten stars from) as M thresh 
=10 x (Mp 0p x /mix)/(/pop x A'), where M pop is the total mass of the Galactic 
population (thin or thick disc), N is the number of stars in the survey and / pop 
is the fraction of stars in the survey that belong to the population. The / m i x is 
the fraction of star forming mass of the population that lies within the survey 
volume. If there is no mixing / mix = M pop /M pop , with M pop being the physical 
mass of the population enclosed within the survey volume. If mixing is maximal, 
i.e., stars born anywhere within the Galaxy can lie in the survey volume, then 
/mix = 1. For the thick disc, M pop =3.9 x 10 9 M 0 and / pop =0.236, which means 
that Mthresh=4.2 x 10 4 M 0 for N=10 6 and we assume / m i x =0.25. Clusters with 
initial masses below 4.2 x 10 4 are outside the detection limits of the survey as 
noted by the green shaded region. Less efficient radial mixing brings stars from a 
smaller number of clusters into our survey volume, moving all of these thresholds 
toward lower mass and making cluster identification easier. 

In Fig. [4j the red dots show the cluster masses from which we would expect 
to recover 20 stars (points on the left) and 40 stars (points on the right) in a 
million-star survey, and the number of such clusters we would expect to find if 
Xmax is 2 xlO 5 (blue curve) or lxlO 6 (green curve). The red error bars show 2a 
Poisson uncertainty on the number and on the size of the recovered groups. We 
see that a smaller survey size would mean fewer stars per formation site, from a 
similar number of formation sites, and severely limit the range of cluster masses 
over which we can explore the ICMF. 

The above calculations assume that clustering exists in chemical abundance 
space and that the clusters are well separated and observational errors are small 
enough (<0.1 dex) such that they can be detected by clustering algorithms. In 
reality, the detectability of clusters in C-space depends upon the dimensionality 
of this space, the intercluster separation and the observational uncertainties on 
abundance measurements. These questions can only be answered with a large 
enough data set that has been homogeneously analysed, such as the GALAH 
survey data. 


6 Discussion 

The argument for a million star survey can be understood most simply as follows. 
Let us assume x max be the maximum mass of a cluster of some galactic population 
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Fig. 4. The number of clusters recovered from a simulated GALAH survey as a function 
of initial cluster mass. The green, dark blue and light blue shaded regions are the ranges of 
initial cluster mass accessible to surveys with 10 6 , 5 x 10 5 and 2.5 x 10 5 stars, respectively. 
The red dots show the cluster masses from which we expect to recover 20 stars (points 
on left) and 40 stars (points on right) in a million-star survey, and the number of such 
clusters we would expect to find if x mELX is 2 XlO 5 (blue curve) or lxlO 6 (green curve). 
The red error bars show the ±2er uncertainty on the number and size of the recovered 
groups. The top curves are for a starburst; the bottom curves are for quiescent star 
formation. 

of mass M. The size of this cluster in the GALAH survey will then be 


— ^max 


/pop-^survey 
/mix fl />op 


( 6 . 1 ) 


For / m i x = 0.1, /pop = 1-0 and M pop = 2.5 x 1O 1 O M 0 , i.e., considering the whole 
galaxy, we have 


n = 


, ^max 

1O 4 M 0 


N* 


survey 


io 6 


( 6 . 2 ) 


So we need N survey to be large so as to observe at least groups of size 10 or more. 

Finally, two factors make it easier to detect groups in thick disc. First, the 7 
is probably higher for the thick disc (-1.5) compared to the thin thin disc (< —2). 
Secondly, the maximum intrinsic size of groups x max is probably also higher for 
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thick disc. Both these effects imply that we will have more large groups for thick 
disc. Additionally, the ratio f P opN survey / M pop is also higher for thick disc which 
also helps. This ratio is purely determined by the selection function and geometry 
of the survey. 

The question remains as to how well we can effectively isolate thick disk 
stars in the GAL AH survey. We will attempt a selection based on a/Fe but 
also explore kinematic selection, and their combination. A differential compar¬ 
ison of the C-spaces for the thick and thin disk datasets may be the most di¬ 
rect route to confirming that higher levels of clumping exist in the thick disk 
(Bland-Hawthorn et al. 2010b) 


Our model is somewhat independent of how the disk formed in that we do not 
need to specify whether the star formation was external or internal to the Galaxyj^] 
In principle, some part of the thick disk may have formed through in situ turbulent 
processes (Lehnert et al. 2014) If the thick disk was formed through accretion, 
it may be possible to detect a flattened dark-matter component, although this 
would be difficult to separate from the baryonic component and the dark halo. 
The process of accretion may indeed have flung open (and globular) clusters into 
the halo and the bulge (Kruijssen et al. 2012)] a signature we can look for in the 
next few years using both Gaia kinematics and chemical tagging. These will be 
chemically distinct from star clusters formed in low-mass dwarf galaxies since most 
dwarfs have mean metallicities well below [Fe/H]=-1. 
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2 For future reference, ‘exogenous’ and ‘endogenous’ may be useful adjectives for describing a 
process, e.g. star formation, that is external or internal to a system. 








